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RO DUCT LON 


Milewciatysis sor Mrcrostrip 1S of great importance as 
Pnis type of mess ion line has found wide use due to 
mesecOompatibility with microwave integrated circuitry. The 
Open microstrip transmission line has been predominant be- 
Gause it can be etched or deposited easily on the substrate. 

Macwcaniyvedialytie work on macrostrip lines has been 
based on a proposed TEM mode of propagation which is, in 
mersense, a Static approximation to a dynamic system. For 
mmeerepentiy tow freauencies the quasi-static theory can 
meecmployed to obtain the characteristics of microstrip 
lines. When the wavelength in a microstrip line becomes 
memparable to transverse dimensions of the line the devia- 
tion from quasi-static behavior becomes significant and 
meet order modes of propagation become possible. 

In this present work a Spectral Domain transform method 
mma pl1eda for Calculating the frequency dependent charac- 
Memrstics of single and coupled microstrip transmission 
mimes. tCitects of geometry on the dispersion and charac- 
Peristic impedance have been analyzed. 

Peounes sh ana 2 Show the geometry of the microstrip 
mene. tlhe strip conductor is assumed to be infinitely 
thin and perfectly conducting and the ground plane and 
dielectric or ferrite substrate are assumed to be infinite 


im extent. 
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Figure 1. Three Dimensional Microstrip Transmission la SNe 


eat 
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Recon 2 


eens cCOuUpLed Microstrip and Regions. 


AZ 





Pew iscmassumea initially that the substrate material us 
lossless and isotropic. Lossy anisotropic substrates are 
mreated later using perturbation theory. 

The Spectral Domain transform method was suggested by 
mean and Mittra, [Ref. 6] and yields a solution to the 
Boundary value problem of microstrip via the method of 
moments. The method was used by Itoh and Mittra to find 
Bie edadispersion characteristics of a single strip. In this 
momk the method has been extended to cover dispersion and 
Baaracteristic impedance of both single and coupled micro- 
mrrmps. itithe casé Of a Microstrip on a ferrite substrate 


mas also been solved using perturbation theory. 
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Ree torekolOn CHARACTERISTICS ON DIELECTRIC SUBSTRATE 


A. FIELD AND BOUNDARY CONDITIONS 
Let the following electric and magnetic fields exist and 
Meepceate in the z-direction as justified by the Hertzian 


Mector Potential functions described in Appendix A. 


o Cee. ¥Z 
E, ko¢ 6 (1) 
Sen .y 2 
Hinman (ac (2) 
Miere Y 1s the Propagation constant. For the lossless sub- 


Strate material 
mee Be 


immometme field expressions in equation (1) and (2) all other 
memponents of electric and magnetic field can be derived 
from Maxwell's curl equations, as shown in Appendix B and 


leads to the following equations. 








ES ie - jou +s ene (3) 
ie coe jwe | el? (4) 
Eas: 5 + jou a a (S) 
Hy = hae; jwe “S| ei (6) 


1A 





Applying boundary conditions at the interface between 
Hme210n 2 and the ground plane, tangential electric fields 


must be zero, so it follows that at y = 0, 


—_ 


iI 
© 


eye 


(7) 


Ey, (6052) 


| 
© 


(8) 


PaIJomaemenceimberazace between regions | and 2, tangen- 
meal electric fields must be continuous: 


At y = 0D, 


E. (x,D,z) (9) 
7 


ee ze) 
ay 


eee) Pere (10) 


Peromenescloetrie fields will exist only in the dielec- 


melee part of the interface and can be expressed as, 


0 Chest rip 
Ee, Ue Ze ee 
al 
ewe elsewhere 
0 Ses ap 
Ee (6,0,2) = C129) 
a 
: Z 
e (xje’ elsewhere 


Poivariy, temgential magnetic fields must be discon- 


menuous by corresponding surface current densities 


ceeaos on strip 


ee sz) = (13) 


—_ 


H. ihe. z) --H 
i. 


0 elsewhere. 


a5 





yz 
a9 , on Serip 


et . = = (14) 


0 elsewhere. 


aa 


Dipset meing tue field expressions of equation (1) 
marough (6) into the boundary condition expressions of 


equations (7) through (14), one obtains 


2 Cc 
me, ?2‘*:9) aoe (15) 
265 863 
Y ax (x, 0) - Te, ay (x, 0) =50 (16) 
k2 95 (x,D) = k?2_ 6° (x,D) (17) 
Cc 1 3 c 3 
1 2 
e h e 
9, (x,D) 9, (x,D) 9, (x,D) 
ax a ay Saar 
7 (18) 
94 (x,D) 
oe? oy 
0 Cles trip 
Be 93 (#0) = (3) 
@ (x) elsewhere 
3 h 0 Chess trip 
94 (x,D) ; 9, (x,D) 
ie ox ee 2H 5, cay 
e, (x) elsewhere 
JX) Ones trip 
k2 gT(x,D) - k? 68 (x,D) = 
Cy Ee? Cy (i (2p) 


0 elsewhere 
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+ jue, —— 


: 964 (x,D) 995 (x,D)| 
oa 


dX dy 
(22) 
ag9(x,D) 845 (x,D) eo 
miimwoces= 9 27 sy | © 
0 elsewhere 


B. SPECTRAL DOMAIN TRANSFORM 
Mme Gveitraletuncc rons MuUSt Satisfy the following rela- 


tion: 


Vio r IS oe) (23) 


where 


Iemma tice =k TRS. (24) 


One introduces the Fourier transform to the a-domain, as 


Suggested by Itoh and Mittra: 
+c 
5 (a,y) = f Daeeyes dx Aer MP Ove. es (e253) 


Thus, the transform of equation (23) becomes 


Peer | 2) + ke FL fo) «0 (26) 
ig ee x * oy? er x 


zZ 
(Jo)? FL f oCxy) 1 + sor F, OGY) 1 
(27) 


2 = 
+ ke B. [ ¢ ] = 0 
where 


d x,y . , ; 
pL. { S000.Y) ] = ja Fr. [ 6(xy) J (28) 


has been applied. 
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Clearly, one can obtain a-domain representation of 


equation (23), 


-a*O(a,y) + = o(a,y) + ki $(a,y) = 0 (29) 
or 
+ = (a? - k2) (a,y). (30) 


Equation (30} should be analyzed carefully both for 
Region 1 and Region 2. 


For Region 1, (free space propagation constant), 


ara Oe 2 ig, 
¥y = Oy om Cotes Ki G51) 
where 
ee ea 
Ky ~ WV OE, and g = ri 
pene a et >)? [ an yee | (Sz) 
1 i! ¢ 


meis the microstrip wavelength and is related to the free- 


Space wavelength A as 


a ee (S35) 


a 
ve 

ie 
Secu t ror ime equation (33) into (32) one can find 


lower and upper values for Yy 


2 7 2m ,2 2 2 
Ciaeu.: (ey Dat X ae 7 2 a. (a4) 
iinieretore a foetivay sear realeaquantity independent of 


Mme Values of a. 
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Pitterariy Lor Region Z, 
Zi r 
get CH fe, - Chea] os 
Be it is clear that Y5 will be imaginary for 


2 2 A 
Tver,” (qryi<a< eo 


and will be real when a is in the range 


ee C57) 
A om A 
wos Bylo) Ve GED ace 


ome coorom i ameme t1ecld functvon differential equation 


molution is as follows: 


Z 
ger 9, (0,y) = yie(a,y) = (a?-KE 8 (ay) (38) 
(s? - y?)$,(s) = 0. (39) 


ioe Lore , 


®,(a,y) = R, (aye V1Y eee: le (40) 


early, tor infinite positive values of y, the field 


Should vanish, or, 


ian o,(a,y) = 0 (41) 
yr 
then, 
(9.7) = Ro(a)e 1” = A(aje “1” (42) 
where 


IMS 





A(a) = R, (a). 


Pet iniveror Region 2a chere are two solutions, cor- 
responding to the real or imaginary character of Yo: The 


solution for Y5 imaginary is 
h h 
o,(a,y) = Bla)Siny,y + C(a)Cosy,y (43) 
where 
(a 
and for Y> real 11s 


o,(a,y) = B(a)Sinhy,y i C(a)Coshy.y. (44) 


f 


Therefore, the transforms of O71. oT O55 $ are given 


by 
Region 1 
ao) 
er(a,y) = A°(aje * (45) 
ot(a.y) = A (aje (46) 
Region 2 
£ Me 
By (o)Sinhy,y + C,, (a) Coshy 5y , Y, Real (47) 
65 (a,y) | 
ie : . : e (48) 
JB. (a)Siny sy + C7. (a)Cosysy » Yo Imaginary 
ep isaiahy vet eogonGartiy Real wie.) 
H 2 H a2 
h 
o5(a,y) 


ae h : _ (50) 
JBp(a)Sinysy + Cia) Cosy5y » 7 Imaginary 
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mimee sper ce1 ipe (e) denotes the electric field case and 


(h) magnetic field case and where 
fifee es ese 
Y9 J¥5: 
specural Domaime representation of the boundary condition 


Sepressions can be obtained by taking the Fourier transform 


Se the equations (7) through (14) with the following results: 


2 e a - 
Betts?) = 0 (51) 
. sah 
-jay?,(a,0) - JOU. By (on 0 je 0 (S22) 7 
2 e _ 2 e / 
ve oi (a,D) = ke 5 (a,D) (53) 
1 2 
e 907 (aD) e 
Joyey(s,D) - juny —gy—— = -jove5 (aD) 
99 (a,D) (54) v 
_ 1W te eee 
Ws sy 
2 @f 2 | 
ig 1 (%?) ~ EB, (a) iy, 
. nah 
Be toD) Jeni ay (%.D) = 2, (a) (56) v 
2 h 7 2 h . 
Be 91640) Kon Og (aD) = J, (a) (s) 
° h Pe 907 (a,D) : h 
ee | ay | (AY P2(e,D) 
. 265 (a ,D) | eal 
ake 5 aan = J_(a). 


hieneomtme derivative transitorm pair 
3 
p { $20ay) ) = -ja0(a,y) 
Mas been applied. 
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Paso tein tnheerield expressions for real (hyperbol- 
ic case) and the imaginary (trigonometric case) values of 
Y> mec tle transtormed equations (51) through (58), the 


following equations are obtained. 


Hyperbolic Case: 


2 e = 4 
ke Cy(a) = 0 . (59) 
2 
: e h : , 
-j LayCy (a) = WH ¥>By (a) J ml (60) 
k? A®(a) = k2 [C&(a)Coshy,D + Bo(o)Sinhy.,D] fae 
F eH 1D H 2 
e h e 
minayA (a) = OUZY4A (ojo: == Jay [By (o)Sinhy,D + 
(62) 
h h 
Cy, (a) Coshy ,D] - juu,y,[B,(a)Coshy,D + Cy(a)Sinhy,D] 
2 € St 
* (a) = EB, (a) (63) 
-j(ayA® (a) - wH,y,A"(a)) = #,(a) (64) / 


5D + Cj (a) CoshyD] = J, (a) (65) 


2 h rn 2 h : : 
ko nen) ko [Bi Co) Sinhy 
1 Z 
5) (0 AP (a) + we,y,A°(a)) + ja feces) Samah Det Ee (ajicosh D] 
J Lay v1 IOV LP Y> H Y> 
-jwe,y,[B,,(a)Coshy,D i C,, (a) Sinhy,D] ee) (66) V 


Similarly for imaginary values of Y5 Dyictkhane ties Sub 
stitutions 
Oda 
Sinhy ,D = jsiny,D (67) 


Za 





Coshy,D = Cosy 5D 


Mmrecduations (59) through (66) a corresponding set for the 
fmeeconometric "case can be obtained. 
Meith con-staies tor clectric and magnetic field expres- 


sions are determined as follows: 





A° (a) = pe F(a) (68) 
Cc 
iL 
oO oo) F(a) (69) 
JWU,Y X WU4Y7K zZ 
ee evel: Cy 

Cr (a) = 0 (70) 
h 

Bp yl) = 0 (7a) 

e _ 1 
Brto) = FR Sinyyp 7260) ey 

LZ 

e Z i 

By (0) = g2-Sinhy,D #2 6%) ey 
Z 
aay = = 1 ea ee ic) (74) 
H WH, Y>,Sinhy 5D Kk? Zs J"x 
Z 
ee | EY (a) - 52, (0) (75) 
It WH, Y'55iny5D 2 Zz Joy 
=p) 


The subscript T denotes the trigonometric case and H 
Memotes the Hyperbolic case. Substituting equations (70) 
mmrough (75) into equation (65) and (66) one can obtain two 


Mees Of equations of the form 


ZS 





F,(a,8)£,(a) + F,(o,8) F(a) = J, (a) (76) 


F,(a,8)E (a) + F,(a,8) Bj (a) = J, (a) (77) 


where 


acer 


memeethe hyperbolic case the a ievewunee Lorn 








k2 k? 
(a, 8) a hy ,D (78) 
F Oh Cotanhy 
ly WUY 4 WY 5 Z 
B aB 
Gee) = | ee + Cotanhy.D (79) 
oH ey OU DY 2 ; 
Pe (a,B) alt PF (a, 8) (80) 
Hi H 
Coeanin as) 
_., | ae oo 2 
(2,8) = - J os 9 Co 
H Pe sc Cc c 
1 ii 2 
(31) 
Gibjiae 
oe Ye 


A similar set fimcdMaerons 10 the trigonometric case 
memebe obtained by use of equation (67). Solutions of 
equations (76) and (77) LO EB, (a) and BL (a) in terms of 
J, (a) and J. (a) may be obtained as follows: 


F,(a,8)d, (a) i F,(a,8)d, (a) 


Ey (a) = 5 (82) 


N 
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-F,(a,8)7,(a) + Fy(a,8)J, (a) 


a (83) 
where . 
Dy = F,(a,8)F4 (058) - Fy(a,6)F5 (0,8), (84) 


Define the following constants, 


F, (058) 
G, (a,8) = Pe cosy 
~F,(a,8) 
G,(a,8) - a (86) 
-F,(0,8) 
G,(a,8) = Lo (87) 
P1083 8) 
pees?) = a (88) 
finally obtaining 
B(a) = G,(a,8)7,(a) + Gy (a,8)7, (a) (89) 
Bo(a) = G,(a,8)J,(a) + 64(a,8)F, (a). (90) 


meee PHYSICAL PARAMETERS 

General equations (89) and (90) show no dependence on 
mee actual physical configuration other than the boundary 
“Conditions. 

Ibpeceneral, the transformed surface current densities 
J, (0) and J, (a) Can be expanded in a set of known basis 


Meaactions as, 





tH 
Ht 8 


J, (a) . (a) (91) 


I) 
lit 8 
Sy 


J, (a) occas (92) 


Since the two conductors are sufficiently:‘narrow, one 
Mmeeeassume that the surface current density in the x-direc- 
tion is zero. By modifying equations (89) and (90), also 
substituting (91) and (92) one obtains 
oo 


leak Dee 


G. (a, 8) (a) 


2 (a) (93) 


Gylo,B) 32 PidaG 


(a) E, (a). (94) 


mireenermore, although one could perform a multiterm ap- 
fameotination to the surface current density in the z-direc- 
faonm, in amie case a one term approximation was found 
Mumpeeralecnt daiter different distributions were tried. For 
meemumnerical results of this present work uniform surface 


current density is used, thus 


Go (a,8)b,j,; (a) = B,(a) (95) 
Gy(a,8)b,j,; (a) = F(a), (96) 

Memcoubstituting equation (92) for the case where i = 1, 
G,(a,8)7, (a) =) (oy) 
G,(a,B)J_ (a) = E (a). (98) 
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It is clear that one only needs to work with either one of 


the above equations. Assuming a uniform current distribution 


= Ti/wW Smee es ta tp 
J, (x) = 


0 elsewhere 


where W is the width of the strip and I, is the magnitude 
@eecne Current, one can obtain the a-Domain Fourier trans- 
omeM , J, (a). iiisewas dene tor three diftrerent configura- 
tions, namely single strip, coupled strips even mode and 
Geupled strips odd mode. 


Romine Seometry im Figure 3 one obtains the following: 


aimgile Strip: 


i W/Z — Sit) 
J, (a) = an e ax J, (a) ~ T. (99) 
OL 
-W/2 a 
Boupled Strips, Even Mode: 
I oy (S/2+W) . 
Ji(a) = a eI "ax + ‘| et %*ax 
-(S/2+W) S77 
(nO oF 
; aW 
J (a) = 2I Cos SEL! es 
Z O 2 W 
(=) 
2 
Odd Mode: 
I =0/2 . -(5/2+W) . 
J, (a) = 7 eFax - | eJ ** Ax (101) 
- (S/2+W) 5/2 


La 








“~~ 


Figure 3. Single and Coupled Strips Even and Odd Longitudinal 
| Mapheii Ui sStrlputions. 
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Sin (ax) 


a(S ai] 
: cv) 


J (a) s j2I Sin (101) 
It is convenient to make equation (98) independent of E (a). 
ies can be done by taking an inner product with a function 
orthogonal to E(a). ADeaiemenis concept, equation (938) 


becomes 


<G,(0,8)7, (a) ,Wla)> = <E, (a) ,W(a)> (202) 


where the inner product is defined by 


++ CO 


Brey Wa)> = J Y(a)W(a) da. (103) 


A suitable weighting function W(a) is the complex con- 
mu@eate of J (a); lee J* (a) or J(-a). 

Then by Parseval's theorem, the right-hand-side of 
eeuation (98) becomes zero because of the orthogonality of 
the E, (a) and J, (a). 


ihereiore equation (98) takes the final form 


+00 
J G,(a,B)IF(a)d, (ajda = 0 (104) 
or 
+00 
/ Gy(a,A/rA")] J, (a) | *da = 0 (105) 
where 


B = 2n/A'. 


a 





pommone cal see the integrand 1s dependent on the ratio 


mummy. OY ertective dielectric constant, om , defined by 
-_ ere 


E.. = ( <t ines 
eff 

Veine the preceding current distributions. for single and 
momopied lines, the phase velocity characteristics can be 
determined by finding those values of the propagation con- 
meome that give a zéro value for the integral at a particular 
frequency. 

Pisce loamethatweonce caquati1on (105) is integrated the 
m@epenaence on the variable a disappears, so one can state 
Mere the whole process is really a function of frequency 
Pememecne Structure’s physical characteristics, namely, width 
Gaetme strips W, and in the case of coupled strips, separa- 
tion S between the Siwiice eilcxmess OL Substrate D, and 
Miemarelilectric'’s relative permittivity e.. 

By observing the behavior of each term equation (105) 
@iencan conclude integration of the IJ, (a) |? term always 
fee a positive area. In the three configuration of 
iimerostrip which were studied, the Fourier transforms of 
pie current distributions always provided a multiplicative 
mietor of (Sino/oa)* which gives rise to a curve having a 


large main lobe and low side lobes in the a-domain. 





Several plots of G,(a,A/A') and the G,(a,A/A")> J, (a) |* 
memm between limits of a; -3000<a<3000 for different A/A' 


waiues 1S shown in Figure 4 and Figure 5 respectively. 
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Figure 4, Gy (a,A/A") versus a. 
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. Figure S. Gyla,rA/A')[F, (a) |? versus a. 
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The family of curves is symmetrical to the y-axis and 
for decreasing values of A/XA', moves along the y-axis with- 
@mtelosing the symmetry. The total area under it monotoni- 
mally decreases. 

The value of A/A' which gives a zero total area under 
maewcuUrve iS the root of equation (105). This gives the 


Sirective dielectric constant e.. and the phase velocity. 
eff 


fee CHARACTERISTIC IMPEDANCE IN TERMS OF THE DISPERSION 
PHARAGTERISTICS 


Mmimcecieraien Ole Can Express the average power in terms 


Smeecne Strip current as, 





ee 
aye 2 aoe c (UG) 
irom which 
2P 
AVE 
tg = TT (107) 


Z 
provided Ls Homdeunigue GCurrent. 
In this study, it was assumed in the dispersion charac- 
memdstic Part that the surface current across the strips 


were uniform; furthermore, the current can be expressed, as 


I -o “dQ = a> ax | (168) 
oO, me, oT 2 


I, a IJ, (x) | dx. (109) 
z JSTRIP 


Mmeene case of a single strip, 


ol 
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J, (x) 


0 elsewhere. 


nor Simplicity ry Pr atoteraraiyveset equal to unity. Thus, 
2 —_— 
De I. (110) 


As stated by Collin [Ref. l], a general expression for 
Eme time average power flow in terms of the electric and 


Meenetic fields is 


— os 
Paypl%Y) = Re i Ir x Hea_da (i) 


where 


E x H*.a 


II 
(Td 


H*. - E mle. 
Z ey, Vie 


This equation becomes 
= y x 
PAVE (x,y) Rel %& J (bao IB H® ) dxdy (eZ) 


Meterrimesback to equations (1) through (6), it is found 
meres the equations for the electric and magnetic fields can 
Memin terms of Hertzian potential Anerson 

poestibstituting equations (3) through (6) ip tne eXpres= 


Sion for average power, 


Payp (Xs) = 74 Re fj (582 7 jupe ed > Dy e58%(. Pas ag 
A (115) 


ae jue2d Je “IBZ axdy + [fae sy + june ie 
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ie: 96° . = j 
(- jas - jue“S— je JB2 axdy eis 


Applying Parseval's theorem 








+00 +00 
h 
. Oo 
PAVE = 7 Re { | abo-(a,y) - jute 229) | 
# e 
-528fo.2) - aweo® (a,¥)] - [58°57 sawneh (ay) 
: (114) 
ie 0° “Ca ) 
Ba ® (oey)) = i ame dady ; 





where the Fourier transform of 3¢/0x has been applied from 


amarron {114). Performing the operations indicated, 
+00 +00 L 
eee 2 & 2 90% (oy) 2 
mun ag -a*Bwe|P(a,y)]|* - w6u| ay | 


e 
- wpe [22 Cea)? - g2uugl eh (ay) |? 
(eles) 
ee es h 
+ jak? (+ ca, 928 Glos) ee) + 9°" (ay) PP My) ; 
. e 302" (a ) h* 90° (a ) 
- jaBb*|o (Coy) ae oe (oy aye | aaay 
where 
eee oe 
) foter ONtaining the général expression for average power, 
Baemedn apply it to both regions for further specific 


equations. 


oe 


ee Oy RS 


® 


ee Re OE 





In Region 1, 


early = 0) 
#(a,y) = A°(aje + (116) 
mee 2) 
a,y) = AMCaje 7 (117) 
965 (a,y) -¥, (y-D) 
—— es eee : (118) 
907 (asy) er -¥,0"-D) 
oT = - yah are (119) 


A 


Substitttineg the preceding equations in (115) one can 


obtain 
+00 +00 
1 “27, (y-D) e 3 
P = —— Re e -a*uwe, |A~ (a) | 
lave At 1 
- Oo D L : 
(120) 


h 
-w8u v2 [AT (a) |? - wBe,¥} [AX (a) |? - aw, BIA COs 


~jak3 (2y, AM (a)Ae™ (a)) + jaB?(2y,A°(a)AP” (a))] dady}. 


Pewucwelear troneete above Expression that one can inte- 
mace With respect to y to obtain 


es a2 + y2 
| 1 1 e 
Pt = OT / a eI Ca) | * 


AVE ee (121) 


+ Bou, JA (a) |? + 20(B2+k2) Re [8° ceoa"*ca} |. 
feces tron 2, clearly there are two expressions for aver- 
age power. When Yo Loeceaeeecduatitonms (4/) and (49) apply 


mmestated earlier. 
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¢ (a) and ph (Gimaue Zero so, 
Teel 


TH 
6©(a,y) = Bo (a)Sinhy eee 
ote y) = C2(a)Coshy (123) 

h 
9, (a,y) h 
a3 = -¥20y (4) Sinhy y (124) 
9O5 (a,y) , 
—sy = -¥2By (4) Coshy py (125) 
and 
D SSS 
eee - 49 2 e 2 
Dy = a, Re If Sinh Yoyla Bwe, |B (a) | 
H HYP” 70 
REGION 


h : h* 
* Bungy51Cy(a)|* + joB?¥,B°(a)Cy (a) 
me GE: (ey idady + Cosh2 

Oo ee y Y2Y (126) 
h | e 
[a*Bwu,| Cp, (a) | * i Bwe, 75 |B (a) | * 
* x 

2 jaB?y,Bo(a)C) (a) - jak3y,C)' (a) Be (a) |] dady 


The y-dependence of the average power in Region 2 dis- 
appears through corresponding integrations as indicated 
below: 

D 


| Sinh2y,D - 2y5D 
Sinh*y,ydy ae 


Roy 
rT (127) 


0 


oy 





D Sinh2y,D + 2y,D 


i] Cosh*y,ydy = ay (128) 
2 
0 
Thus, 
a 
rp = - Re (CSmnhi2yee) ~=5 27. DB) 
“AVE, 167 2 Z 
HYP 
‘ REGION 
a*Bwe 
eZ Cc 2 h 2 
=, IBr(a)|* + Bwu,v,/C, (a) | 


+ 


* 
jaB?Be(a)CAM(a) - jak2C)(a)By,°(a)| do 


, ae ee 
FD 


REGION 


+ 


e a 
Seema tr 30a) |* * joBeB. (a)C,, (a) 


h x 
- Hels, C,, (apiBe (a)ida\ . (129) 


Meets also clear that for imaginary Y5 the average power 
PepressSion can easily be obtained with substitution of (48) 
mao) and the corresponding trigonometric case electric 
and magnetic field constants i.e. equation (72) and (74). 
Peeeetailed exposure to such analysis is given in Appendix 
Cc. 


Then for imaginary Y5 Region 2: 
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1 a* Bwe. 
P. = = Ton Re (2y5D = Sin2y5D) | =a 
AVE. Y> 
TRIG 
REGION 
1! h eC eh 
acon ee ea) |- * ab Bo (ajc. (a) 
Ae Ee caida: + (2y"D + Sin2y"'D) 
Zz.) i 2 2 
TRIG 
REGION 
o* BWH h 2 "! e 2 oa *h 
[ ye Cp Co) | + Bwe 75 | Be (a) | - a8 Be. (a) Co. rc) 
_ 2 h e* 
ak5C. (a) Be (a)ieda } . (130) 


etal power in Region 2 becomes: 


P = Pp +7 GL) 


As strip width decreases one observes a widening a-domain 
meer distribution in Region 1. Also for higher values of 
G@melectric constant where the energy is more concentrated 
[eeor the Substrate, a wider a-domain power distribution is 
@etained in Region 2. 


finally, characteristic impedance can be found as 


a (132) 


mueatlecases, feometry 1s taken into account in the 
@m@erent distribution expressions. 
Plots of a-domain power distribution for Regions 1 and 


eore Shown in Figure 6 and Figure 7 respectively. By 


oo 
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Region 1 Average Power Distribution in o-Domain. 
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Figure 7. - Region 2 Average Power Distribution in a-Domain. 
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marone the inverse Fourier transforms of Py and P, one 
AVE AVE 
could obtain the x-domain power distribution. 


Beno munporeantmtonobserve and prove the continuation of 
Ps and P, curves for the critical points where Y> 
AVE AVE 
Switches from real to imaginary character or vice versa. 
mie enlarged views of Ps and P, given an Figures 8 
AVE AVE 
and 9 clearly show the smoothness of transition from one 


meolon to another. 
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3 Taye 
eo 
ip 
Qo 
Q 
“015 ——(ié‘“‘étC LOCOS 040 z 
xD 





fem: 2. 006 -Oi UNITS INCH 
Meee gers 1) Or +04 UNITS TNH. 
Reg tne 95. P, VieteciSmaU eat) JGhz . 
AVE 
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mime PERIURBATION ANALYSIS OF MICROSTRIP ON FERRITE SUBSTRATE 


Considered here is an extension of the spectral domain 
technique to the analysis of microstrip on a ferrite substrate. 
A familiar perturbation ee ssion hei eie moO pam at VOnmGO ll. 
stant of a hollow, closed boundary waveguide with a change in 


Meterial is given by Helszajn [Ref. 5S] as: 
jw _[ftestaxg) EE + UL [Ax] -H'H*) da 


i) a ee 
{[@ x H'+ ExH*)-a da Glos) 
5 


which primes denote perturbed quantities. This expression 
meemrs Valid for opén boundary structures which support 
bound waves [Ref. 7] as is tne case here. 

Mensider the case of a microstrip on ferrite with mag- 
emeezation perpendicular to the substrate. In this case, 
feeine the lossless dielectric substrate problem is solved, 
there is no change in electric susceptibility so [4x] Or 


The change in magnetic susceptibility is given by 


Xxx 0 Xxz 
[Xn] = joo. 0 (134) 
Xox 0 Koz 


where magnetic loss is permitted. The denominator of (133) 


SO,One Obtamns 


I Fee - beada 
> —— 
oo 7 el 


tT a: [* = JOU, —_—_— nn ‘ (135) 
: oe 


ws 


1s just 4P vc 


45 





Evoldetion of the numerator of (135) leads to: 


+00 D 
*& * 

= ' (x, HH* + x, HH*)dydx 
T'o + T* = joy 

Z 2 ‘foto, (136) 

mere unperturbed fields have been used as a first approxi- 
Meeron to perturbed fields and where because of the geometry 

no demagnetization of r.f. fields takes place. No cross 


terms involving products of H and ie appear in the numera- 


tor of (136) because this product is odd and 


+ 0O + 0Oo 


: flohiogd x= [ Pied ee 0) (137) 
ZX eZ 


= © 
- © 


Referring back to equations (2) and (5), substituting 


mmecdidation (156) ana applying tne Parseval's thecrem one 





Setains: 
: +o if) 
r!o+ re = o = £56 (-jaee(a,y) 
Ll 2 Tt ea ay 
O oO, ane 0 
905 (a,y) 9o° (a,y) 
2 ay ee 2 
mmm en 87 8t.y) = JYE> a 
poe h* ; 
: Xzzkc, 02 ary) % (a,y) > dyda (138) 
where 
Xz ~ Xxx? 


tits proparation constant expression holds for + z 
directed traveling waves depending upon the choice of upper 


M@eerower Sign in (138). 
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Siem wematcempt eto SOlye this equation by defining 
regions of a where Y9 has imaginary and real values which 
are indicated by equations (36) and (37). 

Also substituting previously calculated corresponding 
equations of 67 (a,y) and (35 (a,y))/ay for hyperbolic and 
trigonometric cases into (138) and acca ite Ole em 


dependence the following is obtained: 


r' + fT* = wb Om ee 
167 ie - 
oO ate Y2 


REGION 
h Bes nis 
{ Xacc (-82 a4 + ke cpa)? * jue, Bay! 


° In * ° h* tt 
ee Ge Gm) Ba(e)) - (ue,75)* |B, (a) |? 


—— 1 JOU. Sinh2y.D - 2y50 
167 dal 2 Y5 
hye 


Z 
REGION 


h * 
{ Xaccf(-BPa7 +k e )1Cy Ca) | P25we Bay, (Cy (a) By (a) 


h 


- Ci" (a) BE (a) ) - (we yy)? 135 (0) 172} da. (1e55)) 


Evaluating the difference between perturbed and inper- 
turbed propagation constants, 


r' + T* = Aa + jAB 


(140) 
meee = Ot B' = 4B 


where 
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Noe ead sand AB = 6' - B. 


Pmso, aS stated in Reference 5 


eel 


Xxx = Xxx 7 IXxx- She) 


Substituting (140) and (141) into equation (139) and 
meeenrating real and imaginary parts gives expressions for 
Meond §', loss and phase shift with ferrite respectively. 


Define a common expression 





ct ae JOU, 2y,+sin2y,D 
6 Ti t0r8, 5 
TRIG 
. REGION 
i ees tee |S ae ee 
| Be ; Sinh2y,D - 2y,D 
Z (we ,75) [Ba Ca) | da + Baeette Y 
HYP 
REGION 
h aah e 
(Bhat eke J1Cy Cod 1° + 2we,Bay,(-jC, (a) By (a) ) 
- (we,y,)?| By (a) [7] da (142) 
Then 
Aa = KTyn (143) 
AB’ = Kx!) (144) 


where Mo and Me are defined in Reference 5. 
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boscminedbeper wavelength versus frequency and B/8B 


versus frequency are plotted in Figures 10 and 11. 


49 





996 


lo, = 16.8 


Ul3 


Wi 3.36 


We = o.(00 


QI2 


Osi 


! _— —- a ee 
0Gc UZS g10 O19 g2c 925 


£ Gua) 
PoOerPUNITS INCH: 
F-QO1 UNITS INGH. 


Peeure 10. Loss per wavelength in db versus frequency (GHz). 


50 





oz > 


Z 007 006 QG3 


Uo 


me oOUALES 
fer iL. f=: 


Bacure 11. 


Wyo \6.3 
W = 3-3¢ 





Svemene UNM S INCH. 
eevee UN | > aah. 


B-/8 versus frequency (GHz). 


ol 


agg 


Peur 





IV. COMPUTER PROGRAMMING 


A computer program is written in FORTRAN IV language 
@eweeriyst calculates the phase velocity or A/A' ratio for 
a given geometry of microstrip and parameters of substrate. 
Giemecne program proceeds and calculates the characteristic 
iimeedance for the same set of data using results found in 
the first phase of the calculation namely A/A'. If the sub- 
Beate is ferrite, also phase shift and loss per wavelength 
meeweatculated and curves are generated. 

The program consists of one main routine and three sub- 
fememies. ihe main routine starts by calculating all con- 
eens and parameters wnich are needed te evaluate the 
Various integrals. First equation (105) is integrated and 
Seemuated for an arbitrary initial A/A' value. This process 
Sememues for different values of A/A’ until zero is ob- 
maeed as the value of integral. This computer time consum- 
Mieeebteration process 1s shortened by a simple root finding 
[mene that determines the correct region of A/A' and elin- 
inates others. On the average, after 9 iterations equation 
(105) gives a value which will yield zero value for the in- 
Meera! with an accuracy of 0.01%. The exact number of steps 
will change depending on desired accuracy and physical 
parameters. 

Evaluation of the integral is done by standard IBM 


360/67 system library subroutine DQG24 which uses a 24 


a 





point Gaussian Quadrature formula which integrates poly- 
momaals up to 47th degree. 

This routine has been found to accurately integrate the 
type of curves encountered in this study. 

The DQG24 routine requires an external function sub- 
program which defines the polynomial to be integrated for 
imaginary and real values of Yo: Two different subprograms 
GZR, GZIM were used for this purpose and to test and set 
Maem imits of integration, for real and imaginary character 
of Y meSpectamre ly . 

It is also proven that increasing the limits of inte- 
meatetOn Or using a 32 point integration routine which 
Eveauates polynomials up to 64th degree does not change 
the result more than 0.1%. 

Another subroutine called TRAN supplies the expression 
fOr the transformed current distribution depending upon the 
pouysical configuration of the microstrip. 

As far as the theoretical development and programming 
@emeoncernead the only approximations used are the repre- 
sentation of the axial surface current density by a square 
pulse and the assumption that the transverse surface cur- 
rent density is zero. Different forms of current distribu- 
fons were tried such as 

coshx on the strip 


J(x) = 
et) elsewhere. 
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It was concluded that due to the many zero crossings of 
this transformed current distributven curve and the greater 
number of side iobes which adversely affected programming 
efficiency its use was not justified even though it more 
Mecurately approximates the true current distribution. 

Characteristic impedance was calculated as a next step 
Ppeevaluating equations (129), (130) and (132) using dif- 
Mmaeeiltc Parts OL the same routines used in the first part 
Oeecne calculation. 

PemEnemsUbstrate 15 ferrite the program proceeds to 
another section which evaluates equations (142) through 
mere) and calculates loss and phase shift values as a func- 
muon or Operation frequency. 

meemay De pertinent to quote the typical computational 
mime using the method outlined in this study. 

fecmerne GO step the average time for calculating A/A', 
Zo» and loss and phase shift (if necessary) for a single 
frequency was between 1.5 to 1.9 306 on an IBM 360/67. 

ime computer program accepts the following data: 


= Substrate parameters, ¢e u 


r,” ry - 
- Thickness of substrate, D, in millimeters. 


meiiewseparation between conductors, S, in millimeters 
fonly for coupled strip). 


meimitial frequency (GHz). 

Piewpregram starts calculating values for the initial 
feeauency and repeats the procedure by increasing the fre- 
meemey in 1 GHz steps up to 50 times. Means are also pro- 


vided to compare the half wavelength, A/2, with the thickness 
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of the substrate, D. For 


Se) 


yo 


Bveeprogram Stops due to the presence of higher order modes. 
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VeeOiencont REoULTS ANI) COMPARISONS 


The dispersion characteristic of the microstrip as a 
miferton Of frequency has been evaluated for a variety of 
geometries and these results are presented both in eae, 
and plotted form. In addition, characteristic impedance as 
Meet rOon Of frequency 1S presented for the same geometries 
and dielectric constant. 

meoures 12 through 15 show the frequency variation of 
wemvelength and characteristic impedance of a single strip 
memevarious Values of W/D. Similar families of curves for 
@emptled strips appear in Figures 16 through 19. In this 
case the parameter is S/D and W/D is held constant. 

The results indicate that higher frequencies and wider 
memes increase the effective dielectric constant, ee 
and the characteristic impedance ae of both single and 
weed Strips. The spacing S between coupled strips also 
ieoman effect on ae Increasing spacing, oS. i1nercases 
cr forwoc@eand even modes. The effect on characteristic 
impedance depends upon the mode of propagation. Increasing 
Spacing, S, increases Zo for the odd coupled mode but causes 
a decrease in Zo Lo@uunlemevens coupled mode. | 

All available published data is based on quasi-static 
mimayses Except in a few cases. Basically, in this study 
mmemrra dimension, frequency or D/A, is introduced. Two 


SeepLots are given in Figures 20 and 21 to illustrate this 


Meme. ithe x-axis is D/A, the y axis WePeeand the z-axis 
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meeure 21. Three dimensional plot Lips Teed) wesingle strip. 





A'/A and Zo respectively. One can find the A'/A versus W/D 
and Z versus W/D surfaces already plotted in different 
Mererences, i.e. [Ref. 10] but only for zero frequency. 

This newly added dimension makes it possible to find 
data for a specific frequency in addition to a certain 
geometry. 

In Table I the values found for characteristic impedance 
With two different methods in Reference 3 and generated from 
mais Study are compared. This comparison reflects the agree- 
ment in the low frequency range. 

Also some of the tabulated data in Reference 11] is com- 
[med With data obtained from this study in Table II for 
@oeupled strips. | 

One of the few studies of frequency dependence was by 
Haddad [Ref. 9]. Table III shows data taken from Haddad's 
work and from this study which clearly indicates the agree- 
ment False POuevOweercaucemey ranges. Although Haddad’"s 
melysis was frequency dependent, the data presented by him 
Mes restricted to relatively low frequencies. 

Finally three sample computer output pages are given 
in Tables IV, V and VI. 

Although Pee iawecmnotimand Chagacteristic impedance of 
MeeerOstrips on dielectric substrate depend only upon D/\, 
maemmormalized frequency, this is not true for ferrite sub- 
metees. oilnce the computer program will calculate results 
Mereecither case it is necessary to provide actual dimensions, 


fess > CCC. aS inputs. Where both frequency and D/id are 
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TABLE I. Characteristic Impedances of Microstrip Transmission 
fnessomele Strip. 


“r, Cm 
W/D Z* (ohm) | Z? (ohm) 2° (ohm) 
0.1 135.455 134.352 137.64 
a 2 iso 72 112.255 lec 
0.4 OTT 2 89.909 93.25 
0.7 73.613 71.995 73.53 
120 C2 ais 60.970 64.45 
2x0) 43.149 AGS 100 44.15 
4.0 eee nl 26.027 Dials 
10.0 13.341 12.485 12.58 
oe 16.0 
0.1 85.9659 87.762 oe 
2 71.6954 13025 73.43 
ee 57.4999 58.110 59.20 
7 46.2344 46.217 47.83 
1.0 OS 2 38.948 40.72 
2.0 26.7555 26.248 27.73 
4.0 To? 210 16.300 16.96 
10.0 8.0385 7.807 a 


| * Characteristic impedance obtained by method of moment [Ref. 3] 
[eenaracteristic impedance obtained by conformal mapping [Ref. 3] 


S Characteristic impedance obtained by spectral domain 
method in this study. 
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table Lit. 
SINGLE 
W/D 
ex 
veff 
1GHZ 6.78 
2GHZ 7.05 
ae 
O 
1GHZ 48.52 


Frequency Dependence of Single and Coupled 


Microstrip. 


lew 


et 
Teff 


6.75 


ff -, © ot aw ew 


wmmaddad results 


COUPLED 
(EVEN MODE) 


W/D = 


e* 
veff 


i aoo 


iO Sy ye — 


et 
Teff 


COUPLED 
(ODD MODE) 
We =— t70-S/D 

e* et 
Teff veff 
5.84 5. 85 
So! 5.94 
os Pa 

O O 
Soy 28 39.312 
56. 2Q 39,180 


ieeopectral Domain Transform method used in this study. 
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— 
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maple IV. 


W 


= Imm 


FREQUENCY 


1.000 


2-000 


3.000 


4.000 


5-000 


6.000 


1.000 


8-000 


9.000 


10.000 


11.000 


12.000 


13.000 


14.000 


15.000 


16.000 


17.000 


18.000 


19.000 


20.000 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


D 


= 2mm 


~ D/LAMBDA 


0.00667 


0.01333 


G.02000 


0.02667 


Oe0S233 


0.04000 


0.04667 


“7 a aA Arm 
0.05535 


0.06000 


0.06667 


0.07333 


0.08000 


0.08667 


0.09333 


0.10600 


0.10667 


0.11333 


0.12000 


0.12667 


Velo so5 


2 
2 


L/LPR 


el 


Zetooil 


2.78140 


2.80726 


2283508 


2.86373 


2289254 


2.92108 


2294507 


2-97628 


3.00258 


3-02787 


3205206 


3.07512 


3.09702 


36 lager > 


Solos 3 


Samoa f 


Sauls 1 1 


3618959 


3420464 


12 


W/D = 0. 


LPR/L 


0.36236 
0.35953 
0.35622 
0.35272 
0.34515 
0.34572 
0.34234 
0253909 
Cie33a9S 
0.33305 
Ceooeat 
Cae cic. 
0.32519 
Oe oezes 
0.32074 
0.31874 
0.31688 
0.31515 
0.31354 


0.31205 


5 


Sanprestomputer Output for Single Strip. 


CH IMPEDANCE 


C259 


62./0 


63.07 


Geri Z 


64.68 


65476 


€7.56 


69.51 


1.81 


74.247 


77.51 


80.94 


84.77 


8S.01 


93.66 


$8.73 


104.24 


110.17 


Dre 2 


PZ 1 


OHM 


OHM 


QHM 


OHM 


OHM 


OHM 


OHM 


CHM 


QHM 


QHM 


QHM 


OHM 


OHM 


QHM 


OHM 


QHM 


OHM 


‘OHM 


OHM 


QHM 





Table V. Sample Computer Output for Coupled Strip (Odd Mode) 


W= 1.6mm D = 2.0mm S = 0.8mm W/D = 0.8 5/02 0-4 Ee, = HZ 
2 


FREQUENCY D/LAMBDA L/LPR LPR/L CH IMPEDANCE 
1-000 GHZ 0.00667 2261034 0.38309 39.05 OHM 
2-000 GHZ 0.01333 2.61858 0.38189 38.9% OHM 
3-000 GHZ 0.02000 2.63141 0.38002 38.89 OHM 
4.000 GHZ 0.02667 2.64815 0.37762 38-84 QHM 
5-000 GHZ 0.03333 2.006821 0.37478 38.87 QHM 
6-000 GHZ 0.04000 2<Cotiz 0.37159 39.00 OHM 
~.000 GHZ 0.04667 2.71640 0.36812 39%co OFM 
€.000 GHZ Oe0ss23 2614366 0.3044E 33.609 UnM 
G.000 GHZ . 0.06000 208124S 0.36069 40.31 QHM 

10.000 GHZ 0.06667 2.80248 0635603 41.16 QHM 

11.090 GHZ 0.00253 2.83325 0.35295 42-26 OHM 

12-000 GHZ 0.08000 2.86441 0.34911 43.65 OHM 

13.000 GHZ 0.08667 2.89561 0.34535 45.34 OHM 

14.000 GHZ 0.09333 2.92645 0.34171 47.37 OHM 

15.000 GHZ 0.10000 2-95675 0.33821 49.275 OHM 

16.000 GHZ 0.10667 2.98612 0.33488 5265 2 On 

17.000 GHZ 0.11333 3.01442 0.33174 65.67 Ghit 

18.000 GHZ 0.12000 3.04146 0.32879 59.23 OHM 

19.000 GHZ 0.12667 3206/14 0.32604 63.22 OHM 

20-000 GHZ Wero333 2.0913 9 0.32346 67.65 GHM 





Table VI. 


W=1.6mm D=2.0mm S=0.8mm 


FREQUENCY 


1.000 


2-000 


3-000 


4.000 


2-000 


6.000 


7.000 


8-000 


9.000 


10.000 


11.000 


12.000 


13.000 


14.000 


15.000 


16.000 


17.000 


18.000 


19.000 


20.000 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


GHZ 


D/LAMBDA 


0.00667 
0.01333 
0.02000 
0.02667 
0.03333 
0.04000 
0.04667 
O20 53:33 
0.06000 
0.06667 
0.07333 
0.08000 
0.08667 
0.09333 
0.10000 
0.10667 
0.11233 
C.12000 
0.12667 


Grol33 33 


W/D = 0 


L/LPR 


(Ge 


pecer02 
2.96430 
3.00492 
3.04470 
3.08186 
SelSg oO 
3.14609 
2.17216 
aeons 
gee leas 
3.23740 
3.25422 
3.26921 
3.28263 
3.29467 
3.30552 
3.31533 
See2a23 
3.33234 


3.33975 


70 


S/D = 0. 


== 


0.34164 
0.33735 
0.33279 
0.32844 
0.3244E 
0.32056 


0.31786 


0.31278 
0.31070 
0.30889 
0.30729 
0.30588 
0.30463 
0.30352 
0.30252 
0.30163 
0.30082 
0.30005 


0.29942 


64.22 


64.38 


64.92 


65.82 


67.00 


68.42 


10.00 


41.79 


13.49 


(2435 


77.26 


19.220 


81.16 


83.214 


85.13 


Gee 3 


89.14 


Sie 5 


S307 


95.18 


Sample Computer Output for Coupled Strip (Even Mode) 


OHM 


JHM 


QHM 


OHM 


OHM 


OHM 


OHM 


CHM 


OHM 


OHM 


CHM 


QHM 


OHM 


CHM 


OHM 


OHM 


CHM 





shown as printed output the calculation is based on the sub- 
strate thickness, D, provided as an input. A width of D = 


2mm has been used in the calculations presented here. 
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VI. CONCLUSIONS 


In this thesis a hybrid mode analysis of single and 
coupled microstrips (of equal width) on dielectric substrate 
Daiemdescribed. The resulting equations were solved by the 
method of moments applied in the spectral or (Fourier) trans- 
form domain. A computer program was developed to calculate 
mevelength and characteristic impedance versus frequency. 

The original results obtained by this method were compared 
mietme low frequency limit with the quasi-static results 
Ot several other investigators and were found to be in ex- 
eerrent agreement. 

mie method was extended to microstrip on ferrite where 
the propagation constant was calculated using perturbation 
mmeory . Nurtemetes1 CVWORvtehOmwOheresultinge €Xpressions 1S 
again carried out in the spectral domain. A computer program 


has been developed for this analysis also. 
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APPENDIX A 


PUR Y VECTOR POTENTIAL FUNCTIONS 


Pee MAGNETIC HERTZIAN VECTOR POTENTIAL FUNCTION 
fonsctaecr a homogencous, source-free, isotropic region; 


Mmemee, there 1s no charge density and 


V°eE = 0 (Al) 


Maem vector Algebra, it 1S known that the divergence of 
@yecLtor is zero if the vector is, in turn, the curl of 


another vector; therefore, one can state, 


VeVxT, = 0 OT, (A2) 


— 


| ~jwuvxmy . Cy 
For time-varying fields 
VxH = jwe E (A4) 
momeapplying equation (A3), 


VxH = jwe (-jwuVxm, ) 


=e Ven 


h 
ee. bas (AS) 
= k Vx (TT, + k Vo) 
H = kom, + Vo (A6) 
mere , 
Ie = wpe. 
Pamilarly, for time-varying fields, 
VxE = -jwpH (A7) 
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SO, applying equation (A3), 


Vx (-jwuVxT, ) = -jou(k?m, + Vo) 


- V*7, = k*T,_ + Vo. 


h h 


| 
<] 


a 


(A8 ) 


(AQ) 


Up to this point $ has been defined; arbitrarily choose 


d = Vem: 


iaen, equation (A9) becomes, 
i Demer g2) 

V Tat k i 0. 

Also, it is known that 

V°B = uVeH = 0 

uVe (k*m, + Vo) = 0 

Kem, + VeVo = 0. 


Summarizing, for TM modes, 


E = -jwuVxm, 

+7 _ 125 =.= 

H Kem, + VVermy 
= VXVXT) 


semidarly, for TE modes, 


H = jweVxT, 
E = k*n7_ + VWVet 
e 
= VxVxT.. 
e 


B. TE AND TM MODES FROM VECTOR POTENTIALS 
1. TE Modes 


E = -jwyvxm 
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(A10) 


(A11) 


(A12) 


(A13) 


(A14) 


(A15) 


(A16) 


(A17) 





y= kén + VVet 
= Vxvxnt 
where 7 Saelstiles: 


y2qe ~ Kae = 0 


For TE modes, E. = 0; therefore 
ame an. 
me iat 5 - sy 


from where, 


ie 
is. = Ty, = 0 
Let: 


7 = pe Y2a 


Z 
Dien. £rom@equation (A13) 
h 2,h 
V76 + k?26° = 0 
t? e? 
where, 
Pee YK? 
Summarizing, 


a k29"(x,y)e*¥? 


H, = tye V2y gn 


t 
Be = a a xH 


tr3| 
| 

a 

NO 
= 
+ 
<] 
4 


o 
tt 
de 
Ee 
m 
<] 
* 
~ 
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(A18) 


(A19) 


(A20) 


(A21) 


(A22) 


(A23) 


(A24) 


(A25) 
(A26) 


(A27) 


(A28) 


(A29) 





For TM Modes, He = 0; therefore 


T° = o°(x,yje "7a, (A30) 
and, - 

Vig" + k29° = 0 : | (A31) 
where, 

ke = y2 + k2, GRSZ) 


Summarizing, 


e tvzZ 
E = k?¢ ese (A33) 
= —- ertyz 
E, = tyV,¢°e if (A34) 
— = = 1WEe a — 
ie + — aE. (A35) 


oommonse@eey and) (2) corresnond to equations (A33) 


ome (AZ5), respectively. 


vo 








APPENDIX B 


TRANS ERSESELECTRIC AND MAGNETIC FIELDS 


a 


Consider Maxwell's point form of Amere's law, 


VxH = 


ee 
ay 


0H 
OZ 


rad 


oH 
al. ° 


OX 


jweE 
dH 7 
OZ Jee 
= jwe(E,a 
3 
dH 
—~ = JweE, 
JH, 
9x = se 
dH 
Bay = Jwek, 


and, due to the 


can apply 


@ 
Lo 


lL 
a 
' 


(B1) 


(B2) 


(B3) 


(B4) 


(BS) 


PacopomeemecrOrative f1elds, 1.¢,, ee one 


JweE 


JWEE 
ey 


Jwek. 
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(BO) 


(B7) 


(B8) 


(B9) 





Similarly, starting from Maxwell's point form of Lenz' 


law, 
VxE = - jw (B10) 
JE JE JE JE JE JE 
= Ze y = x Z =: va 4 x 
pr dy Oz ) Ee oy ( 02 dX es ao ( dX oy 
eecrira? Hay ala) Gem? 


irom where, 


9E 

ay yE, = -jwuH, (B12) 
aE. 

ee = Jour (B13) 

JE aE. 

sk - S = -jontt, (B14) 


Supstituting the value of ES. from equation (B13) into 


equation (B7), one obtains 


oH pee OE. 
Eh ee a —— - ywyH 
gy YAY = yO ax 7 Jeu 
3H 9E 
24,2 = ie ee 

H (y?tw?ue) = y 2 - jue 2 (B15) 

Hee OE 
My ~ Kz OY ~ay ~ IPE OX 


Similarly, substituting the values of EY from equation 
mee) into equation (B8), of Hee from equation (B8) into 
Semation (B12) and of Hy Proneecauiactone(h/) ICO Equation 
woes), one obtains 


- For H., 
x 
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3H 1 oL. 
ie | ox | (CME YC oy ee 
oH JE 
2 2 re 1 eee 
LY +4) eo Y —xx + Jue Ty 
1 3H, dE. 
H. = Re ( yx * jwe iy ). 
- For E_, 
jE. 1 JH, 
me yy Coax * Jeeey) 
JE 30H 
zZ 2 ae 1 eee 
EY Cy Ge en) ey xy * JMU x 
1 JE. 3H 
Hy ~ Kz ay * JME ax). 
et Or E> 
; JE, oH. 
fe Y 5 = WJ ay JweE ) 
SE, 3H 
Ee CY Boe e) = ¥ as Jwu “Dy 
1 or. 3H 
Ee Kz TORS - JjJWU Oy Nee 
Summarizing, 
1 dL. oH, 
a = K2 CGas = jou sy ) 
1 0b. 3H 
By ~ Kz CY ~gy * JOH ox 
1 oH, 3E 
an Ke Sema oy 
1 oH ob. 
Hy “ Kz OY “gy 7 JRE 3x 
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(B16) 


(Cem) 


(Sere) 


(B19) 


(B20) 
(B21) 


(B22) 





Equations (B3) through (B6) correspond to equations 


woo) through (B22), respectively. 
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APPENDIX C 
AVERAGE POWER EXPRESSION IN REGION 2 


(TRIGONOMETRIC CASE) 


Electric and Magnetic field potential functions and 


derivatives with respect to y become 


®5(a,y) = jBp(a)Sinyyy (C1) 
h ? 
O7 (ay) = Cpa) Cosyyy (C2) 
905 (a,y) 
= ae iY 5Bp (a) Cosyyy (C3) 
h 
9. (a,y) 
ay = 3 GC) Sinygy. (C4) 


Bupstituting anto (115) 


D 
e : " 
P, ae Re I T2806 1850) Siny3y 1-08 
th TRIG ~0 
REGION 


| - cr 7(a)Siny! Paes 31 5B. (a)Cosy! ale WU 581 C2 (a) Cosy) alee 
+jak3[ (Cr (a) Cosysty) (- jyyBo (a)Cosy5y) + (- jBC (a) Siny! oY) 
(-¥3Cq(a)Sinyyy)] - ja8?[(JB2(a) Sinyty) (-y¥Ch” (a) Siny#ty) 

e (Cc r " (a) Cosy" ¥) GY5Br (a) Cosysy)]|dyda (CS) 


Then, 
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D 
P, = zh Re If Sin? yydy[-a?Bwe, | By (a) |? 
TRIG~ JO 


REGION 


-wBuy? [Cr pla) |*-a°k5y5 ich Na) Be (a) -a8? Y5 ‘Be S(a)cn (a) Jda 


D 
J Cos*yyydy[- whe, ys? |By (a) | 2-27 wuy8| Cr (a) | 
TRIG 


REGION 


+ ak5y' Honk Na) Bo oe: ap? yoBy. (a) cn” (a)}]da). feo) 


Also the integration with respect to y in Region 2 may be 


accomplished analytically as indicated below: 


2 2y4D - Sin2y¥D 
Sin*ytydy = —+——,-—— UGE) 
Z 4y5 
0 
: 2y4D + Sin2yyD 
Cos*ytydy = —-——_,--—_— (C8) 
Z 4y5 
0 


By substituting (C7} and (C8) in (C6) and rearranging 


one can easily obtain 


. 1 Bwe 
P, ae aa Re [ (2y, Des Sinz y’ a a 2 [Be (0, ) |? 
TRIG 


REGION 


+ BWH,Y>5 c (a) | 4ak2Cr (a) Be” (cat aB2 BO S(a)Cr (a) J da + 


alba 
(275) + Sin27! etna [Cy ey eo 57) |B (c))| - 
TRIG "2 


REGION th h k 
Bee Gonie. (a) - is Ge IB (a) ]da). (C9) 


oD 





- The final equation (C9) corresponds to equation (130). 
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